Debye-Onsager-Relaxation-Effect beyond linear Response and Antiscreening 

in Plasma Systems 

K. Morawetz 
Fachbereich Physik der Universitdt Rostock 
18051 Rostock, Germany 

The quantum kinetic equation for charged particles in strong electric fields is derived 
and analyzed with respect to the particle flux. It is found that the applied electric field is 
screened nonlinearly. The relaxation field is calculated completely quantum mechanically 
and up to any order in the applied field. The classical limit is given in analytical form. 
In the range of weak fields the deformation of the screening cloud is responsible for the 
Debye-Onsager relaxation effect. The result beyond linear response presented here allows 
to investigate a field regime where no screening cloud is present. The descreening field 
is determined as a function of thermal energy density of the plasma. For stronger fields 
the moving charge is accelerated by accumulated opposite charges in front of the particle. 
This can be understood in analogue to the accoustic Doppler effect. A critical field 
strength is presented up to which value a thermalized plasma is possible. The range of 
applicability of the treatment is discussed with respect to applied field strength and space 
gradients. 



I. INTRODUCTION 



High field transport has become a topic of current interest in various fields of physics. In semiconductors 
the nonlinear transport effects are accessible due to femtosecond laser pulses and shrink devices [jp. In 
plasma physics these field effects can be studied within such short pulse periods One observable of 
interest is the current or the electrical conductivity which gives access to properties of dense nonideal 
plasmas |jj . In high energy physics the transport in strong electric fields is of interest due to pair creation 

In order to describe these field effects one can start conveniently from kinetic theory. Within this 
approach the crucial question is to derive appropriate kinetic equations which include field effects beyond 
linear response. The most promising theoretical tool is the Green function technique Jl],|5|-[7j . The resulting 
equations show some typical deviations from the ordinary Boltzmann equation: (i) A collision broadening 
which consists in a smearing out of the elementary energy conservation of scattering. This is necessary 
to ensure global energy conservation (ii) The intracollisional field effect, which gives additional 
retardation effects in the momentum of the distribution functions. This comes mainly from the gauge 
invariance. 

One of the most important question is the range of applicability of these kinetic equations. Up to which 
field strengths are such modifications important and appropriate described within one-particle equations 
? In j|] this question has been investigated for semiconductor transport. It was found that for low fields 
ordinary transport equations are sufficient while for high fields the field effects are absent. The critical 
value for this range is given by a characteristic time scale of field effects tJ, = mh/(eE ■ q) compared with 
the inverse collision frequency. This criterion is a pure quantum one. It remains the question whether 
there are also criteria in the classical limit. For a plasma system we will find that there is indeed a critical 
value of the field strength which can be given by classical considerations. 

At low fields we expect the linear response regime to be valid. Then the contribution of field effects 
to the conductivity can be condensed into the Debye- Onsager relaxation effect [^0f^5| which was first 
derived within the theory of electrolytes |l6|-|20[] . Debye has given a limiting law of electrical conductivity 
|]l6| which stated that the external electric field E on a charge e is diminished in an electrolyte solution 
by the amount 
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where k is the inverse screening radius of the screening cloud. This law is interpreted as a deceleration 
force which is caused by the deformed screening cloud surrounding the charge e. Later it has been shown 
by Onsager that this result has to be corrected to 

E(l ^= ) (2) 

if the dynamics of ions is considered. 

In this paper we will give the complete result beyond linear response for the static case (|l|) 

E(1-^F{E)). (3) 

We find that F(E) represents a nonlinear dependence of the Debye- Onsager relaxation- effect on the 



field strengths (45), which is a result beyond linear response. We see that the electric field is limited 
below a value of 

E<-. (4) 

e 

Above this limit no quasi- equilibrated transport is possible, i.e. no thermal distributions are pertained in 
the system. Then we have to take into account nonthermal field dependent distributions which have been 
employed to study nonlinear conductivity pl"|Jj2| . In particular new types of instabilities are possible 

We shall give in this paper a regime below this limit where the surprising effect of antiscreening can 
occur. We shall see that in the range 

kT kT , . 

— > E > 0.7 — 5 

e e 

the sign of F(E) in (||) can change which means that already the field blows away the screening cloud. 
However, the charged particle feels an accumulation of opposite charges in front of its path. This can be 
considered as an analogue to the Doppler effect. We will call it Charge - Doppler effect (CD- effect). This 
has consequences on the conductivity which has to increase within this range. If the field reaches values 
of E = 0.7 kT '/ 'e we have no screening cloud. This effect of descreening has been described by different 
methods . 

First we give a short rederivation of the field dependent kinetic equations from the Green's function 
technique. In the third section we shall derive the field dependent current analytically. This will produce 
the main result of this paper which is the nonlinear relaxation effect. The complete quantum case is 
given up to cubic response. The classical result is calculated analytically for complete response up to all 
orders in the field strength. In the last section we shall discuss some applications of this CD- effect. 

II. REAL-TIME GREEN'S FUNCTION TECHNIQUE 

We would like to consider a system of charged particles of fermions or bosons under the influence of 
an external constant electric field. Coupling the field in vector potential gauge we have the Hamiltonian 

H = J2 I dr* i *(r > *)£r i (5v r --A(t))* i (r,t) 

J 1 c 

(6) 
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With the Coulomb potential 

^•(r-r) 



e 2 ZiZj 
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A. Definitions and Equation of Motion 



In order to describe correlations in highly noncquilibrium situations, we define various correlation 
functions by different products of creation and annihilation operators 

G>(1,2) = < *(l)*+(2) > 

G < (1,2) = <*+(2)*(l) > . (7) 

Here <> is the average value with the unknown statistical nonequilibrium operator p and 1 denotes the 
cumulative variables (ri, s\, t\...) of space, spin, time etc. 

From the above definitions we can build a causal function with the help of the step function 

G(l, 2) = e(h - t 2 )G> (1,2) - e(t 2 - h)G< (1,2). (8) 

We use the equation of motion of the creation and annihilation operators to derive kinetic equations, 
which may be solved with the appropriate choice of boundary and initial conditions. Applying the 
equation of motion for the field operators in the Heisenberg picture, one finds the famous Martin Schwinger 
equation hierarchy p3), where the one-particle Green's function G couples to the two-particle one G 2 , 
etc. 

A formally closed equation can be reached with the introduction of the self-energy 

J d2V(l-2)G 2 (12,l'2+) = j dlE(l, 1)G(1, 1') . (9) 

c 

In order to obtain solutions and the path of integration G, it is necessary to specify the initial conditions. 
In many physical situations, the condition of weakening of initial correlation is an appropriate choice 

lim G 2 (12l'2')| tl = t ,+ 6 =G(11')G(22')-G(12')G(21') (10) 

t\ — >— OO 

which leads into a vanishing selfenergy in (|^). Therefore we can choose for the path C the ordinary time 
integration along the real axis if we subtract just terms which ensure this limit. The result reads 

+ 00 

J dlX(l,i)G(l,l') = J dl{S(l,l)G(l,l')-E<(l,l)G>(l,l')} 

C -oo 

/ti pti fOO POO 

£>(1,T)G<(I,1')+ / £ < (1,1)G<(1,1') + + / E<(1,1)G>(1,1')- / £<(1,1)G>(1,1') 
-oo J ti J t'^ J — oo 

(11) 

which indeed vanishes for t\ ^ t\ t = — oo. This can be alternatively formulated as a contour of time 
integration known as Keldysh contour. This means with other words that the weakening of the initial 
correlation which breaks the time reversal symmetry is equivalent to the Keldysh contour integration. 
With the expressions (||) and ( |ll|) we can finally write the first equation of the Martin-Schwinger hierarchy 
in the form of the Kadanoff-Baym equation p6| , p7l f7| 

- i(G^G< - G<G - 1 ) - i (G R Z< - £<G A ) - * (S^G< - G<£ A ) (12) 

where the retarded and advanced functions are introduced as A R (1,2) = —iQ(ti — t 2 )[A > ± A K ] and 
^4^(1,2) = i<d(t 2 — ti)[A > ± A K ]. Here operator notation is employed where products are under- 
stood as integrations over intermediate variables (time and space) and the upper/lower sign stands for 
fcrmions/bosons respectively. The Hartree- Fock drift term reads 



G^ll') = ( ih± - - S^(ll') ] 5(1 1') (13) 



with the Hartree Fock self energy 
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E HF (1,1')= (T5(ri-r[) J dr 2 V{r 1 -r 2 )G < (r 2 t[r 2 t 1 ) 

+V(n - ri)G < (r 1 t 1 r' 1 t' 1 )) S(h - t[) (14) 

where G(r 2 ,ti,r 2 ,ti) = n(r 2 ,ti) is the density. 

Introducing Wigner coordinates T = (t\ + t 2 )/2,r = t\ — t 2 ,etc, we finally obtain the time diagonal 
part as [E8| 



d_ 
~df* 



OO 

fw(p,T) = J dr [{G > b,T-^-r),E<(p,T-^r)} - {g<( P , T - I -r), E>(p, T - I r)} 



(15) 



Here fw(p,R,T) — ^fig < (p, R,T,t — 0) denotes the Wigner distribution function and {,} is the anti- 
commutator. This equation is exact in time, but according to the assumed slowly varying space depen- 
dence we have used gradient expansion for sp ace ya riables and dropped all R-dependence for simplicity. 
This criterion is discussed in the last section ( [II C| ). 



B. Gauge invariance 

In order to get an unambiguous way of constructing approximations we have to formulate our theory in 
gauge invariant way. This can be done following a procedure known from field theory f29fl . This method 
has been applied to high field problems in ]3C[ |. 

One can introduce a gauge-invariant Fourier-transform of the difference coordinates x 



g(k,X) = J dx expi ^ J d\ + -A^{X + Xx)} 



}g{xX). (16) 



For constant electric fields, which will be of interest in the following, one obtains for the generalized 
Fourier-transform 



g(k,X) = j dx e^ x » k " +erET ^g(x,X) 



where the x function was chosen in such a way that the scalar potential is zero = (0, — cET). Therefore, 
we have the following rule in formulating the kinetic theory gauge-invariantly 

1. Fourier transformation of the difference- variable x to canonical momentum p. 

2. Shifting the momentum to kinematical momentum according to p = k — eET . 

3. The gauge invariant functions g are given by 

g(p, T) = g{k - eET, T) = g(k, T) = g(p + eET, T). (17) 

We shall make use of these rules in the following sections. In ]3l|] this procedure is generalized to two- 
particle Greens functions and leads to the field- dependent Bcthe- Salpeter- equation. 

C. Spectral function 

The spectral properties of the system are described by the Dyson equation for the retarded Green 
function. A free particle in a uniform electric field, where the field is represented by a vector potential 
E(t) = — - A(T) leads to the following equation 
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.^-e(p--A(t)) 



(18) 



This equation is easily integrated 



G$(p,tt') = -iO(t-t')exp 



\Jdu e(p - ~i4(u)) 



(19) 



For free particles and parabolic dispersions, the gauge invariant spectral function 



follows 



Ao(k,u)) = 2 / drcos 



2mh 



e 2 E 2 

"i 

24mh 



o 



2?r^. fk 2 /2m-Tiuj 



(20) 



where Ai(x) is the Airy function ]34[ and = (h 2 e 2 E 2 /8m) 1 / 3 . It is instructive to verify that (pcf) 
satisfies the frequency sum rule. This interaction-free but field-dependent retarded Green's function 
can be obtained from the interaction-free and field-free Green's function by a simple Airy transformation 
p5[ . This is an expression of the fact that the solutions of the Schrodinger equation with constant electric 
field are Airy-functions. The retarded functions can therefore be diagonalized within those eigen-solutions 
p6| , p0| . It can be shown that ( p(i| ) remains valid even within a quasiparticle picture |35| ], where we have 
to replace simply the free dispersion k 2 /2m by the quasiparticle energy efc. 



D. The Problem of the ansatz 

In order to close the kinetic equation (filf), it is necessary to know the relation between G > and G < . 
This problem is known as an ansatz and must be constructed consistently with the required approximation 
of self-energy. The conventional way to do this is to change the correlation functions into the generalized 
distribution function and into the spectral one, which is still an exact transformation 

G<(pwi?T) = A( P luRT)F{ P ujRT) 

G>{ P ljRT) = A(puRT)(l =p F{ P loRT)), (21) 



with the spectral function from eq. (|20|). 

Assuming the conventional ansatz, we replace the u> dependence of the distribution function by their 
quasiparticle value, which is given by the quasiparticle distribution function / 

G<(kLuRT) = A(kujRT) f(kRT) 

G > {kujRT) = A(kojRT) (1 t f(kRT)). (22) 

This is quite good as long as the quasi-particle picture holds and no memory effects play any role. As 
we shall see, the formulation of kinetic equations with high fields is basically connected with a careful 
formulation of retardation times. Therefore, the simple ansatz, called KB ansatz fails. 

Another obscure discrepancy is the fact that with the old ansatz, one has some minor differences in 
the resulting collision integrals compared with the results from the density operator technique. With the 
old ansatz, one gets just one half of all retardation times in the various time arguments [ p8|]32] |. This 
annoying discrepancy remained obscure until the work of Lipavsky, et al. []37f where an expression is given 
for the G < function in terms of expansion after various times. We can write in Wigner coordinates 

G<(p,T,t) = f(p,T-^)A(p,T,T). (23) 
The generalized- Kadanoff- Baym (GKB) - ansatz is an exact relation as long as th e self energy is taken in 



Hartree- Fock approximation. Together with the requirement of gauge invariancc (II B) the GKB ansatz 
finally reads 
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G<(k, T, t) = f(k - T - M) A(fc, r, T) 

G>(fc,T,r) = (1 Tf(k- ^,T- M)) A(fc,r,T). 

(24) 

At the end, we can use the spectral function derived in the preceding section (|2Cj) to obtain the resulting 
ansatz valid for high applied electric fields 



G < (krRT) = e -*( ekT+ ^ r3 ) f (h 



eEM \t\\ , N 

-^l,R,T-U-\. (25) 

In order to get more physical insight into this ansatz one transforms into the frequency representation 

oo 

G < (kujRT) = 2 / drcos ( cor - e(k, R, T)- - ^r 3 j f(k-—,T--). (26) 



Neglecting the retardation in / one recovers the ordinary ansatz (|22|) with the spectral function (J20|). 
The generalized ansatz takes into account history by an additional memory. This ansatz is superior 
to the Kadanoff-Baym ansatz in the case of high external fields in several respects p5[ : (i) it has the 
correct spectral properties, (ii) it is gauge invariant, (iii) it preserves causality, (iv) the quantum kinetic 
equations derived with Eq. ( |3l| ) coincide with those obtained with the density matrix technique |28j , and 
(v) it reproduces the Debye-Onsager relaxation effect JTH) . 



E. Quantum kinetic equation in high fields 



With the help of the gauge invariant formulation of Green's function ( |II B| ) , we can write the kinetic 
equation (rL5|) in the following gauge-invariant form 



oo 

A /a(fc; T ) + eEV k f a (k, T)=J2[dr[lG>(k- ^r, r, T - 1), E< (k - ^r, -r, T - 1] 

b o 



G< (k - ^r, r, T - I), E> (k ^r, r, T- T -) 



(27) 



where all functions are now gauge-invariant. This kinetic equation is exact in time convolutions. This 
is necessary because gradient expansions in time are connected with linearization in electric fields and 
consequently fail pq |. The gradient approximation in space has been applied assuming slow varying 
processes in space. This local dependence of all functions on R is suppressed in ( p7j ) and furtheron. Wc 
will discuss the validity in the last section ( III C| ). 

> 

To obtain an explicit form for the kinetic equation we have to determine the selfenergy E < . For the 
selfenergy we use the statically screened Born approximation 

Z< b (pRTT) = ^f^(P + P' - P - P"')K 2 b ( P - P) 

x G> (p', -r, R, T)G< (p, r, R, T)G< (p>, r, R, T) 

(28) 

with the potential given by the static Debye shape 
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and the static screening length k. The latter expression is given by 



c e ° lc 

in equilibrium and nondegenerated limit. Here s c is the degeneracy for specie c and T c the corresponding 
temperature. 

Introducing ( p8| ) into (|27]), the kinetic equation for two charged particle scattering with charge and 
mass rrn in high electric fields reads p5p^ , |35| ] 

d d 

b 

_ 2s b f dk' a dk h dk' h 

t/2/1 ,m [i f/ / /\ T Et 2 /e Q k a e 6 k b e a k a e h k b \l 

x K fc (ka - k a ) / dr cos \ (e a + e b - e a - e b )~ - — + ^ f 

J L h 2h \ m a m b m a m b / ) 

o 

X {/^(l - - fb) - fah(l - f' a ){l - ft)} ■ 

(31) 

Here we have written, e.g.. f b for f b (k b — e b Er,T — r) for simplicity. If we had used the conventional 
Kadanoff and Baym ansatz (g2j) we would have obtained a factor 1/2 in different retardations |2Sj| . This 
would lead to no relaxation effect at all jlq ]. Furthermore it is assumed, that no charge or mass transfer 
will occur during the collision. Otherwise one would obtain an additional term in the cos function 
proportional to r 3 . This kinetic equation is derived in second Born approximation but with complete 
field dependence. Generalizations can be found for the T-matrix |HJ approximation resulting into a field 
dependent Bethe-Salpeter equation or for the RPA approximation J35| resulting into a Lenard-Balescu 
kinetic equation. 

Two modifications of the usual Boltzmann collision integral can be deduced from ([si]) : 

(i) A broadening of the ^-distribution function of the energy conservation and an additional retardation 

in the center-of-mass times of the distribution functions. This is known as collisional broadening 
and is a result of the finite collision duration (3^] . This effect can be observed even if no external 
field is applied. It is interesting to remark that this collisional broadening ensures the conservation 
of the total energy pOfl . If this effect is neglected one obtains the Boltzmann equation for the field 
free case. 

(ii) The electric field modifies the broadened 8- distribution function considerably by a term proportional 

to t 2 . This broadening vanishes for identical charge to mass ratios of colliding particles. At the 
same time the momentum of the distribution function becomes retarded by the electric field. This 
effect is sometimes called intra-collisional-field effect. 



III. FIELD EFFECTS ON CURRENT 

We are now interested in corrections to the particle flux, and therefore obtain from (|l|) the balance 
equation for the momentum 

^ < k a > -ne a E = ^2< kJ a B » > . (32) 

b 

For the Boltzmann collision integral the momentum transfer on the right side would be zero. Here we 
shall find a finite value which will be represented as renormalization of the external field E similar to the 
Debye-Onsager- Relaxation field in the theory of electrolyte transport [T^-[2C|] . This effect can be shown 
to be a result of the deformation of the two-particle correlation function by an applied electric field. 
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To proceed we assume some important restrictions on the distribution functions. First we assume a 
nondegenerate situation, such that the Pauli blocking effects can be neglected. In order to calculate the 
current for a quasistationary plasma we choose Maxwellian distributions 



(33) 



with the thermal wave length Af = 2Trh 2 / (rriiTi) and a partial temperature for species i which can be 
quite different e.g. in a two - component plasma. 

Then the momentum conservation in (31) can be carried out and we get for momentum transfer 



< k/|f >= - 



2s b f dqdQdk 2 ■.[,,,, „ , 



X cos 





Eq 


efc e a 


t 2 




" ~2h 


mb m a _ 





{/„(k)/ 6 (Q + f ) - + q)/ 6 (Q - | ) } 



(34) 



where we have shifted the retardation into the distribution functions. The second part of the distribution 
functions can be transformed into the first one by k + q — > q and Q — q/2 ^ Q + q/2 with the result 



< kl% b >= 



x dr cos 



2s b 



dkdqdQf b (Q)f a (k)V 2 (q)q 
kq qQ 





2{ih m a h mbh , 





Eq 


efc e a 


T 2 




" 2h 


mb m a _ 





(35) 



with the reduced mass fi 1 = l/m a + l/mj,- The angular integrations can be carried out trivially and 
we get 



< pit >= |a 



m b m a 



2 

. q t , 



(36) 



with js(x) = (a; cos a; — sinx)/x 2 . The two integrals over the distribution functions li can be done with 
the result 



h\a\ = 



qr 



dkf a (k) sin(-^-) 
m a n 



S3„ _2„ 2 «a m „ 



(37) 



and correspondingly I 2 [b] . 

We now introduce dimensionless variables 



q = 


2y^ fJ ,T 




2Tt 




t = 








h 




T = 






2 \m a 


4- mb 






e b e a 


e = 


4T 3 / 2 


TOb m a 



-T, 



(38) 



and obtain 



h = 8_n^T J dtis{yt 2 e)sHy 2 t)e -v> 



Using now the screened Debye potential ( p^ ) we finally obtain 

%n a n h (? a e\sb T 

2 l — J 3 



(39) 



7V 2 
1 e o 



h = dz 



(z 2 + l) 2 




j, ■ i ,2\ sm ( 2;2 ^) -z 2 l 2 
dljs(xzl) — \ — -e 



Therein we used y = zb and I = tb with the quantum parameter 



and the classical field parameter 



e E 

X = T 



m a m b 

-eft e, 



b 2Tn \m a + mi, m a + mb 



(40) 



(41) 



(42) 



With this form (40) we have given an extremely uscfull representation because the field effects, contained 
in x, are separated from the quantum effects, which are contained in b. Both expansions will be performed 
and discussed. To proceed we give a series expansion in the field parameter x 



00 7 _|_ 1 

T - 2 Vf-l) H1 — — T 2W f 



r2fc + l 

J 3 



fc=0 



dz- 



2fe+l 
3 



Z 



2fc+4 



(z 2 + l) 2 




(2fc + l)! / , .: 



2-2fc 



dz- 



2 J ~ (z 2 + 1) 
o 



2 lFl (2 + 2fc, -,- — ). 



The last integral can be done analytically and is treated in the next section. 



A. Classical limit 



(43) 



First we give here an exact expression for the classical limit. We obtain from (|43| ) for 6^0 

4 « = (^J^^ = ^ ( _ 1)t(1 _ 2t , 



(44) 



With the help of this result we can sum (|43|) with the restriction x < 1 from (42) and obtain 



'3c 



2k- 1 



8 ^ 2fc + 3 

k=0 



-y 



,2fe+l 



irx , . 
'24 F ^ 



F(x) 



4- 3x 2 



1 — X 2 X 



-artanhx 



(45) 
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This is the main result of our calculation. It gives an exact classical result for the field dependent Debye- 
Onsager- relaxation- effect up to field strengths x < 1. 

Introducing the classical result (El) into (EH) we find from ([36|) and (p2|) the following relaxation field 



at V E 

with 



d ( 5E 

< k a > -n a e a E 1 + — ^ ] = (46) 



8E a _ e a 7T \ - 4rt b efc — 



Sb 2 ^ m \2 F ^)- (47) 



We see that for a plasma consisting of particles with equal charge to mass ratios, no relaxation field 
appears. The link to the known Debye- Onsager relaxation effect can be found if we assume that we have 
a plasma consisting of two species with opposite charge but equal masses and temperatures T a = T& = T. 
Then @) reduces to 

S ^L = ^ F{Xc) (48) 
E 6e T y c> K ' 

with 

eE , v 

x. = 5£ (49) 

and F[x) — l + o(x 2 ) from (^5|). This formula together with the general form ( ^7|) is the main result of the 
paper. It gives the classical relaxation effect up to any field strength and represents a result beyond lincarc 
response. Further it leads to a natural limitation of our treatment with respect to the field strengths 
which will be discussed in section ( [II C| ) . We like here to point out that the opposite special case of light 



particles screened by heavy ions lead to the same relaxation effect, because the limit m{,/m a — > oo leads 
to (EJ|) with — x c . Because of the even character of F(x) this leads to the same expression as in the case 
of identical oppositely charged particles. 

B. Quantum correction 

The complete quantum case of ( |43| ) can be given by performing the integral. The result, which can be 
verified by MATHEMATICA reads 



+ (-l)M2£ + l)! (^3 1 F 1 (2 + 2fc,|^) + 2b 2 (l + k) lFl (3 + 2fc,|^) 

(50) 

With the help of this integral (^FJ) we have given the complete quantum result of the series expansion 
with respect to electric fields It is easy to create any kind of degree of response by this way. Only 

odd exponents of the fields appear in ( |43| ) which gives an even series for F(x) in ([45|). In the following we 
give only the first two parts of the expansion with respect to the field. The quantum effects are included 
completely. 

1. Quantum lineare response 

The first quantum correction reads 
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rln m T (, 1)2 b r- £ b2 \ - A 

/ 3 1 (fc = 0)=g l + T --^e, - + - erfc(-) 



(51) 



Renaming the quantum parameter 6/2 from (^) into x, the result of [ jT^j is reproduced by a different 
way of calculation. Unfortunately in |lE{| a misprint occured where a factor 6/2 in front of the er/c was 
missing. 

2. Cubic response quantum result 
As a next term we can give the cubic Debye- Onsager field effect in complete quantum form 



Jf (jfe = 1) = — — (15V^+ 186i - 1186? - 726? - 8 
dV ; 646i K 

V^re 6 * (-15 + 156? + 1506? + 766? + 86*) erfc(6i)) 



(52) 



with 6i = 6/2. 

In the Fig. g a we plot the quantum versus classical result for linear response and cubic response versus 
the quantum parameter (|4l|). We see that the cubic response is less influenced by quantum effects than 
the lineare response result. The general observation is that the quantum effects lower the classical result 
for the relaxation effect. A detailed analysis of quantum effects on the lineare response can be found in 
II- 




FIG. 1. The ratio of quantum to classical Debye-Onsager relaxation effect (ml) versus quantum parameter b 
of ([ll]). In the left hand figure the lineare and cubic response result, 1$ and /fof ( ji^ ) is plotted separately. In 
the right hand figure we give the relaxation effect up to cubic terms in field parameter x of ( fi^ ) for different field 
strength. 

InFig. |l|b we give the ratio of quantum to classical result for the relaxation effect up to cubic terms 
in fields for different field strengths represented by values of x from (^). We see that the quantum 
effects become more important with increasing field strength. The effect of sign change for field strengths 
x > 0.7 described above can also be seen in the quantum effects at certain values of 6. 



C. Validity of the Calculation 

We have seen that the electric field is limited to values x << 1 for x from For equal masses and 

temperatures of plasma components this condition translates into 
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(53) 



In this limit F(x) has a singularity and no quasi equilibrated transport is possible, i.e. no thermal 
distributions are pertained in the system. Then we have to take into account nonthermal field dependent 
distributions which have been employed to study nonlinear conductivity [^l|^3| . 

The condition (^33j) allows for different physical interpretations. Within the picture of the screening 
cloud we can rewrite (|5^) into 

v 2 

eE < m-&. (54) 
td 

This means that a particle moving on the radius of the screening cloud rjy — 1/n with thermal velocity 
v 2 h = T/m should not be pulled away by the acting field force. This can be considered as a condition for 
the possibility of a screening cloud at all in high fields. For fields higher than this limit we will not have 
screening at all. We can discuss this limit also via the energy density which can be reached in a system 



by the applied field. We can reformulate once more the condition (53) to find equivalently 

e E 2 



4tt 



< nT. (55) 



This means that we have essentially nonthermal effects to be expected if the energy density of the field 
becomes comparable with the thermal energy density. 

The validity criterion (|53| ) can now be used to check the weak space inhomogeneity which has been 
assumed during our calculation. Quasi- equilibrium in a plasma system with fields can only be assumed 
if the field current is accompanied by an equivalent diffusion current 

dn 

jficid = euEn = -jdiff = eD— (56) 

dx 

using the Einstein condition \i = eD /T one gets 

eE = T^. (57) 
n dx 

Combining this elementary consideration with our condition ( |53"|) we obtain a limitation for space gradients 

dn 

< n (58) 



d (ki) 

where our treatment of field effects and kinetic equation is applicable 



IV. DISCUSSION 

It is interesting to discuss the classical nonlinear relaxation effect via the function F(x) in fl45|). In 
figure ^ we have plotted the corresponding dependence. One recognizes that the linear response regime 
is possible until x w 0.2. Then the relaxation effect becomes lower with increasing fields until it leads 
to a complete vanishing at x = 0.69609. At this point no screening cloud is established and therefore 
no relaxation effect present. If we increase the field we see that the effect changes the sign. This range 
means that a charge is accelerated more than the applied field would allow for because of the surrounding 
opposite charges. We call this region antiscreening. It occurs in the range 

^L >E >0.7— . (59) 

e e 

Antiscreening means that the field is already blowing away the screening cloud. However, the charged 
particle is feeling an accumulation of opposite charges in front of its path. This can be considered as an 
analogue to the Doppler effect. We will call it Charge - Doppler effect (CD- effect). This has consequences 
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on the conductivity which has to increase within this range. Of course, the total energy of the system is 
not changed, because this acceleration comes from thermal energy of other particles. 



1.0 



— 0.0 



-1.0 



0.0 




descreening 
0.69609 



1.0 



FIG. 2. The scaled nonlinear Debye-Onsager relaxation effect ( p7| ) or (^) versus the field parameter (|42[). The 
lineare response regime is given by small numbers of x where no remarkable deviation from classical Debye result 
(|l]) is seen. The screening cloud is deformed further with increasing fields until no relaxation occurs at x ~ 0.7. 
This is complete descreening. Further increase of fields change the sign of the relaxation effect and the moving 
particles become accelerated instead of decelerated. This range is indicated by antiscreening. 



The interpretation of accumulated charges in front of a moving charge is underlined by consideration 
of a very heavy projectile with mass m a which goes through a plasma of lighter particles mi,. Then 
we get from ( |47| ) also a sign change. We can imagine this effect as the accumulation of plasma charges 
in front of the moving heavy ion. This, of course, gives a field additional to the external field. As an 
application of this nonlinear Debye Onsager relaxation effect we may think of stopping power of a heavy 
ion in a plasma system |41|. Within the discussed special field range of the CD-cffcct we expect a slight 
acceleration of ions and reduction of the stopping power. This is especially of importance for current 
experiments, where heavy ion beams are stopped in a nonlinear plasma as target. 
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